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NOTATION 


Density  gradient  in  surrounding  fluid,  a  = 


d  dz 


Initial  buoyancy  parameter  (theory).  Equation  [44] 
Vorticity  mixing  coefficient.  Equation  [9] 

Half  distance  between  cores  of  vortex-pair 
Stratification  parameter  (theory).  Equation  [45] 
Constant 

Energy  dissipation  coefficient 


Energy  dissipation  parameter.  Equation  [46] 

Total  energy  of  convected  mass 

Experimental  buoyancy  parameter,  Equation  [44] 

Geometrical  parameter,  J  *  0  for  planar  geometry, 

J  =  1  for  axial  symmetry 

Virtual  potential  energy  coefficient.  Equation  [38] 
Virtual  mass  coefficient 


Virtual  kinetic  energy  coefficient.  Equation  [37] 
Vertical  component  of  total  momentum 

Parameter  defined  by  Equation  l 4?) 

Radial  distance  from  center  of  rising  mass,  r*  » 

Mean  radius  of  rising  mass 

Non-dimensional  mean  radius  of  rising  mass,  R  •  R/Rq 
Experimental  density  stratification  parameter. 
Equation  ( 33) 

Time 

Non-dimensional  time,  t  *  W  t/z 

o  '  o 

Time  at  which  the  maximum  height  of  rise  is  reached 
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u,v,w  Velocity  components,  see  Figure  14 
W  Vertical  velocity  of  rising  mass 

W  Non-dimensional  vertical  velocity  of  rising  mass, 

W  =  W/W 

'  o 

# 

W  Vertical  velocity  of  ideal  vortex-pair 

z  Height  of  rising  mass  center  above  its  virtual  origin 

z  Non-dimensional  height  of  rising  mass  center,  z  =  z/z 

zmax  Maximum  height  reached  by  rising  mass 

0  Modified  entrainment  coefficient 

7  Local  vorticity.  Equation  [7] 

T  Total  circulation  about  a  single  vortex 

£,tj,C  Coordinate  system,  see  Figure  14 

p  Local  density  inside  convected  mass 

Pe  Density  of  surrounding  fluid 

P1  Average  density  of  convected  mass.  Equation  [43] 

Ap  Density  difference,  Ap  *  p^-  pg 

^  Velocity  ratio, 

Subscripts 
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ABSTRACT 

A  theory  for  the  motion  of  two-dimensional  turbulent  vortex 
pairs  in  homogeneous  media  has  been  developed  based  on  separate 
velocity  scaling  of  the  Internal  and  external  flow  fields  in¬ 
volved  in  the  motion  and  taking  into  account  variations  in 
volume,  circulation,  momentum,  and  energy.  Based  on  the  results 
obtained  from  this  theory  (I)  a  simplified  theory  (II)  is  de¬ 
rived  to  deal  with  the  rising  motion  of  turbulent  vertex  pairs 
in  stratified  media.  The  theoretical  results  are  compared  with 
systematic  experimental  observations. 

In  theory  (i)  the  ratio  of  internal  to  external  velocity 
scales,  f,  is  introduced  as  an  Important  variable  and  the  theory 
is  specifically  derived  for  the  two  limiting  cases  of  weak  ($  -*  1) 
and  strong  (y-  »  l)  circulation.  The  weak  circulation  theory 
leads  to  results  similar  to  those  obtained  in  the  past  using 
theory  based  on  complete  similarity  and  momentum  conservation; 
i.e.,  z  ~  t^.  The  strong  circulation  theory  leads  to  results 
which  depend  very  much  on  the  way  in  which  vorticity  from  the 
shear  layer  is  invested  into  the  vortex  pair.  When  ingested  so 
as  to  cause  annihilation  (cancellation)  of  the  ingested  vorticity, 
the  asymptotic  trajectory  ia  z  t^.  Under  the  same  conditions 
the  velocity  ratio,  i>,  increases  toward  an  asymptotic  value,  and 
the  virtual  momentum  coefficient  for  the  motion  tends  to  zero. 

As  a  result,  the  asymptotic  motion  (assuming  vorticity  annihila¬ 
tion)  corresponds  to  a  motion  with  complete  similarity  and  with 
energy  conservation. 
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A  comparison  of  experimental  observations  of  rise  versus 
time  and  radius  versus  height  with  theory  (l)  lend  strong  support 
to  the  strong  circulation  theory  and  suggest  that  ingested  vor- 
ticity  may  be  largely  annihilated. 

Based  on  these  finding  for  homogeneous  flows,  a  simplified 
theory  (II)  for  stratified  media  was  developed  upon  the  assump¬ 
tions:  (i)  the  motion  is  determined  by  conservation  of  volume, 
mass,  and  energy  (neglecting  vorticity  and  momentum);  (ii)  com¬ 
plete  similarity  (dR/dz  =  0,  a  constant).  Good  agreement  was 
found  between  the  predictions  of  this  theory  and  the  results  of 
systematic  experiments,  and  particularly  for  the  maximum  rise  of 
height. 

I.  INTRODUCTION 

Ideal  Vortex-Pairs.  Plow  visualization  studies  carried  out  by 
Scorer  ( 1957 ) j  Woodward  (1959)  and  Richards  (1965)*  indicate  that 
the  shear  layer  which  is  formed  between  a  moving  isolated  mass 
of  fluid  and  the  stationary  surrounding  medium  tends  to  roll  up 
and  create  a  flow  field  which  resembles  (in  two  dimensions)  the 
one  associated  with  two  line  vortices  of  equal  strength  but  op¬ 
posite  sign,  separated  by  a  distance  2b,  so-called  "vortex-pairs." 
The  possibility  of  vortex-pair  motions  in  an  inviscid  fluid  was 
considered  and  analyzed  over  100  years  ago  by  Sir.  W.  Thomson 
(1867).  His  analysis  applies  only  to  an  idealized  vortex-pair 
in  which  each  vortex  has  a  highly  concentrated  core  which  is  set 
into  motion  only  by  the  influence  of  the  other  vortex.  Such  an 
ideal  vortex-pair  moves  through  the  surrounding  fluid  in  a  di¬ 
rection  perpendicular  to  the  plane  joining  the  vortex  cores  and 
with  a  velocity  W*  determined  only  by  the  pair  separation,  2b, 
and  the  circulation  about  a  single  vortex,  T,  according  to  the 
relation 
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A  unique  feature  of  this  idealized  vortex-pair  motion  is 
the  existence  of  a  closed  streamline  and  a  finite  captured  mass, 
as  indicated  by  the  oval  in  Figure  2a.  Thomson  (1867)  calculated 
the  semi -axes  of  the  oval-shaped  captured  mass  to  be  2.09b  and 
1.73b  so  that  the  cross-sectional  area  is  approximately  3.62irb2 
and  the  ratio  of  width  to  thickness  is  1.21. 

Under  certain  circumstances  it  is  entirely  possible  that 
carefully  balanced  vortex-pairs,  approximating  Thomson's  ideali¬ 
zation,  can  be  formed.  The  motion  around  the  vortex  centers  must 
be  affected  by  viscosity  in  a  real  fluid,  but  as  long  as  the  vis¬ 
cous  cores  do  not  extend  close  to  the  bounding  closed  stream¬ 
line,  the  flow  within  and  without  this  streamline  may  be  so 
closely  matched  that  no  large  shearing  motions  or  accompanying 
drag  are  associated  with  the  motion  of  the  captured  mass.  In 
fact,  nearly  ideal  vortex-pairs  are  sometimes  found  in  the  wakes 
of  lifting  surfaces.  Figure  la,  and  are  known  as  "contrails", 
see  Scorer  (1958)  and  Spreiter  and  Sachs  (1951).  Of  course,  the 
concentrated  vorticity  in  the  vortex  cores  tends  to  diffuse,  and 
does  so  rapidly  when  the  flow  in  the  core  is  turbulent. 

Turbulent  Vortex-Pairs.  The  probable  short  lifetime  of  ideal 
vortex-pairs  under  turbulent  conditions  gives  special  importance 
to  vortex-pairs  whose  behavior  is  governed  by  turbulent  entrain¬ 
ment;  indeed,  it  is  these  kinds  of  motions  which  are  most  com¬ 
monly  observed  in  nature,  as  in  the  case  of  a  mass  of  fluid 
forced  out  rapidly  through  an  aperture.  Figure  lb,  or  in  the 
convection  of  isolated  masses  in  nature.  Figure  lc,  or  in  the 
bent-over  and  rising  chimney  plume. 
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Turbulent  vortex-pairs  are  characterized  by  the  fact  that 
the  interior  motion  does  not  match  the  outer  flow  at  the  boundary 
of  the  captured  mass,  so  that  a  region  of  high  shear  exists  there, 
accompanied  by  the  production  of  vorticity  and  by  turbulent  en¬ 
trainment.  In  other  words,  these  vortex-pairs  move  with  a  ve¬ 
locity  W  not  equal  to  the  velocity  W*  derived  from  Thomson's 
model.  Equation  [  1 ] . 

We  may,  in  principle,  generalize  Thomson's  model  to  con¬ 
sider  those  cases  where  the  velocity  of  translation,  W,  has  a 
more  general  relation  to  the  velocity  W*  which  characterizes  the 

internal,  rotational  motions  of  the  vortex-pair.  Two  distinct 

* 

cases  suggest  themselves,  in  theory.  In  one  case,  W  >  W  ,  and 
the  convected  mass  loses  volume  to  the  surrounding  fluid  and 
continually  shrinks  in  size;  we  denote  this  vortex-pair  as 
"underdeveloped",  see  Figure  2b.  In  the  other  case,  W  <  W*,  and 
the  "overdeveloped"  vortex-pair  gains  mass  through  the  entrain¬ 
ment  of  exterior  fluid.  Figure  2c.  Of  these  it  is  the  latter 
motion  which  is  most  commonly  observed  in  nature  and  forms  the 
main  subject  of  this  work. 

Within  an  overdeveloped  motion,  the  velocity  at  the  boundary 
inside  the  vortex-pair,  as  seen  by  an  observer  moving  with  it, 
will  be  larger  than  the  velocity  of  the  surrounding  fluid  Just 
outside  the  boundary  of  the  pair.  Accompanying  the  velocity 
gradients  thus  created  across  this  boundary,  shear  stresses  are 
exerted  by  the  vortex-pair  on  the  surrounding  fluid,  resulting 
in  the  entrainment  of  outer  fluid  and  a  general  increase  in  the 
volume  of  che  convected  mass,  see  Figure  3.  Within  the  high 
shear  zones  at  the  boundary  on  either  side  vorticity  of  sign 
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opposite  to  that  within  the  respective  interior  is  created,  and 
is  pulled  down  around  the  bottom  of  the  rising  mass  toward  the 
plane  of  symmetry.  To  the  extent  that  the  ingested  vorticity 
remains  on  its  own  side  of  this  plane,  the  vorticity  within  the 
interior  will  be  steadily  reduced ;  of  course,  ingested  vorticity 
of  opposite  sign  does  have  a  chance  to  mix  and  thus  to  annul  it¬ 
self,  depending  on  the  efficiency  of  mixing.  Should  effective 
annihilation  of  ingested  vorticity  occur,  then  the  initial  total 
vorticity  within  one  side  of  the  pair  would  be  conserved  in  time. 

As  for  the  kinetic  energy  implicit  in  the  motion  of  the 
vortex  pair,  it  must  be  continuously  reduced  with  time  due  to 
turbulent  dlssi'ation. 

Self-Similarity  in  Vortex-Pair  Motions.  It  is  a  striking  char¬ 
acteristic  of  free  turbulent  flows  in  homogeneous  media  at  suf¬ 
ficiently  high  Reynolds  numbers  that,  under  similar  circumstances, 
the  flows  at  different  points  in  space  or  time  can  usually  be 
reduced  from  one  to  another  upon  normalization  by  an  appropriate 
length  and  velocity  scale  (self-similarity).  This  is  true,  for 
example,  of  the  flow  at  different  downstream  sec  a.ons  of  turbulent 
jets  and  wakes.  It  is  therefore  natural  to  expect  that  a  turbu¬ 
lent  vortex-pair  exhibits  complete  self-similarity  during  its 
life  time,  and  this  assumption  has  been  made  in  all  theoretical 
treatments  of  the  subject,  starting  with  Morton,  Taylor,  and 
Turner  (1956).  Two  important  consequences  of  this  complete  simi¬ 
larity  are:  (l)  conservation  of  the  ratio  of  internal  and  ex¬ 
ternal  velocity  scales,  W/W*,  during  the  motion;  (2)  linearity 
of  the  length  scale  of  the  convected  mass  with  the  distance 
traveled  frdffi  a  virtual  origin. 
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This  latter  result,  predicting  that  the  traces  of  the  side 
boundaries  of  the  convected  mass  form  a  wedge,  is  independent 
of  the  dynamics  of  the  motion  and  serves  to  provide  a  check  on 
self-similarity.  In  fact,  a  number  of  previous  experiments  on 
self-convecting  masses  claim  to  confirm  this  behaviour  to  a  rea¬ 
sonable  approximation,  see,  e.g.,  Scorer  (1958),  Woodward  (1959) 
and  Richards  (1965). 

It  is  obvious  to  ask  whether  a  "natural"  value  of  the  ve¬ 
locity  ratio  W/W*,  or  the  same  thing,  of  the  constant  £  *  dR/dz 
is  observed,  independent  of  the  original  circumstances  giving 
rise  to  the  convected  mass.  The  answer  seems  to  be  no.  In  the 
present  experiments,  two  distinctly  different  ranges  of  value  of 
dR/dz  differing  by  a  factor  2,  have  been  repeatedly  measured; 
these  correspond  to  two  different  stroke  lengths  in  the  ap¬ 
paratus  used  to  originate  the  vortex  motions.  Furthermore,  al¬ 
though  the  present  data  may  be  claimed  to  correspond  "in  a  reason¬ 
able  approximation"  to  a  constant  value  of  dR/dz,  yet  quite 
consistent  deviations  from  linearity  e^ist  between  the  traces  of 
part  radius  and  distance  traveled,  see  Figure  4.  These  deviations 
are  such  that  dR/dz  seems  actually  to  Increase  throughout  the  ob¬ 
served  motions. 

In  view  of  these  facts,  and  for  other  reasons,  it  seems 
desirable  to  attempt  a  more  general  theory  of  the  motion  of  tur¬ 
bulent  vortex  pairs,  based  on  the  assumption  of  separate  velocity 

£ 

scaling  of  the  internal  and  external  motions;  i.e.  allowing  W/W 
to  vary  continuously.  Afterwards,  a  simplified  theory  pertaining 
to  motions  in  stratified  media  will  be  developed  and  the  results 
compared  to  experimental  observations. 


\ 
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II.  THEORY  (HOMOGENEOUS  FLOWS) 

Separate  Similarity  of  Internal  and  External  Flows.  We  visualize 
the  vortex-pair  motion  to  be  divided  into  internal  and  external 
flow  fields,  separated  by  a  thin  region  of  high  shear,  which  also 
forms  the  boundary  of  the  captured  mass,  see  Fig.  3.  We  assume 
that  each  flow  field  is  itself  self-similar. 

Internal . 

w^x^t)  =  Wi(t)-wi  j|  ;  jj  [2] 

External. 

we(x,y, t)  =  We(t)-we 


and  similarly  for  the  other  velocity  components. 

We  let  ip  -  W  /W  ,  where  ip  is,  in  general,  not  constant  in 
time  as  it  is  in  the  case  of  complete  similarity. 

We  choose  W^  as  the  circumferental  velocity  averaged  over 
the  inner  boundary  of  one-half  of  the  vortex  pair  and  Wg  the 
same  except  averaged  over  the  outer  boundary.  (The  inner  and 
outer  boundaries  are  separated  by  a  thin  shear  layer). 

Volume  Changes.  The  volume  of  fluid  comprising  the  vortex  pair 
increases  continuously  with  time  due  to  entrainment  into  it.  Be¬ 
cause  of  the  similarity  assumed,  the  rate  of  entrainment  of 
volume  must  (in  two  dimensions)  be  proportional  to  a  character¬ 
istic  velocity  and  a  characteristic  length.  We  take  for  the 
former,  the  velocity  difference  W^W 

=  2ttR(W1-  We)-a (ip) 


[4] 
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or,  W 

^  (*-l)-a(*)  =  (*-l)a'(f)  [5] 

Note  that  the  dependence  of  the  proportionality  constant  a'  on 
the  velocity  ratio  ^  has  been  left  unspecified. 

Circulation  Changes.  The  circulation  T  about  one  half  of  the 
vortex  pair  is,  on  account  of  the  inner  similarity,  proportional 
to  the  length  scale  and  inner  velocity  scale.  On  account  of  the 
way  in  which  W  was  defined, 

T  «  RW  [6] 

The  fluid  entrained  into  the  vortex  pair  from  the  surround¬ 
ing  high  shear  layer  carries  vorticity  of  opposite  sign  to  that 
already  within  the  interior,  see  Fig.  3  and  thus  reduces  it  to 
the  extent  it  is  not  annihilated  through  mixing  with  vorticity 
being  entrained  in  the  opposite  lobe.  The  strength  of  entrained 
vorticity  must,  on  account  of  similarity,  be  proportional  to  the 
ratio  of  a  pertinent  velocity  and  length  scale.  In  particular: 

(V  V 

Y(entrained)  «  - g -  [7] 

The  flux  of  entrained  vorticity  takes  the  form: 

Vorticity  Flux  «  Volume  Flux  •  y (entrained) 

«  [27rR(Wi-  We)-aty)] 


[8] 
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Finally,  the  change  in  circulation  may  be  related  to  the  flux 
of  vorticity: 

=  -Ai • 2ma(^ ) • (W^-  Wg)2  [9] 

or, 

d(W  R) 

— =  -Ai  a'.(Wr  We).(^-1)  [10] 

The  value  of  the  parameter  Ai  will  depend  in  part  upon  the 
extent  to  which  entrained  vorticity  from  each  side  mixes  to¬ 
gether  causing  annihilation.  In  the  case  of  complete  annihila¬ 
tion,  Ai  s  0. 

Momentum  Conservation.  The  momentum,  M_,  of  the  vortex  pair 

—  Z 

is  conserved  in  motion  through  a  homogeneous  medium.  It  may  on 
account  of  similarity  be  expressed  in  the  form, 

M 

—  =  K  ty)*W  R2  =  const.  [11] 

p  mvr/  e 

where  K(vO  is  a  momentum  coefficient.  Its  form  may  be  deduced 
through  use  of  the  identity.  Lamb  (19^5)*  pg.  229> 

M 

-f-JJy  c<UdC  [12] 

The  latter  integral  can  be  taken  separately  over  the  interior 
and  shear  layer  of  the  vortex  pair,  which  yields  terms  upon 
making  use  of  similarity, 

//  y&W  «  w^2  [13] 

int . 

//  ■yCdedC  «  (W  .  W^R2  [14] 

shear  layer 
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As  a  result,  the  form  ofKj^)  consistent  with  our  assumptions, 
is  seen  to  be, 

Kffl(*)  =  Ki-  K [15] 

where  Ki  and  K2  are  undetermined  constants. 

Energy.  The  total  kinetic  energy  in  the  vortex  pair  motion, 
taking  account  of  similarity,  may  be  expressed  as  the  sum  of  two 
terms, 

=  K,W.2R2  +  K  W  aR*  [16] 

p  I  i  e  e 

while  the  dissipation  takes  the  form, 

d  22 

Si - W13R  •  CD(*>  1171 

where  for  large  values  of  Cp  must  approach  some  limiting 
value  CD',  while  for  small  values  of  *(*-  1).  CD  —  CD(1)- 

Laws  of  Motion.  Limiting  Oases;  Weak  Circulation  (i>  ll-  In 
this  case  the  inner  and  outer  flows  are  almost  matched  and  the 
deviation  from  the  ideal  vortex  motion  is  small.  The  laws  of 
motion  in  their  appropriate  form  become, 

-  (^-l)*  (l)  (5*) 

d(WR) 

— - - -  Ava'(l)-(*-l)s.We  (10ft] 

K  ( 1 ) ' W  R2  »  const, 
m'  '  e 


Volume. 


Vortlcity. 


Momentum. 


[  Ha] 
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Energy. 


K(l) 


d (Wg2R2  ) 
dz 


W 


Combining  [5a]  and  [10a]  leads  to  the  result. 


We  *  )  Al  *  0 


Whereas,  [lla]  requires 


W  «e  —  or  z  «  t^ 
6  R2 


[17a] 


[18] 

[191 


so  that  the  presumed  motion  can  take  place  only  if  Ai  =0  or 
Ax  *  (^-l^"1.  Combining  [17a]  and  [10a]  leads  to  the  require¬ 
ment  that  the  velocity  ratio  be  constant  and  have  the  value 


♦ 


[20] 


At  the  same  time,  dR/dz  is  required  to  be  constant  and  to  have 
the  value. 


*  -  1 


CD(1) 

Sk^TT 


[21] 


We  leave  till  later  a  discussion  of  comparison  with  experiment, 
but  we  may  note  now  that  the  prediction  [19]  is  similar  to  that 
of  the  previous  theory  based  on  complete  similarity. 

Strong  Circulation  (♦  »  l).  In  this  case  the  interior  cir¬ 
culation  is  very  strong  relative  to  the  ideal  value  and  the 
deviation  from  the  ideal  vortex  motion  is  large.  The  appropriate 
laws  of  motion  are, 


Volume. 


dR 

dz 


♦*o' (* ) 


1 5b] 


1 
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Vorticity. 

d(W.R) 

dz  -  -  A1a'(^)-W1-^ 

[10b] 

Momentum. 

(Kx-  Kzi/ ) •  W^R2  -  ^.(Mz/p) 

[lib] 

Energy. 

K«f  (W.R  )*  *=  -  W  *-R-C  ' 
e  dz  i  '  i  D 

[17b] 

Combining  [10b]  and  C 17b 3  leads  to  the  requirement  that  a(^)  be 
constant  and  equal  to, 

C  ' 


2K  -Ai 
e 


[21] 


Since  similarity  requires  that  W &/W  be  constant,  we  may  here¬ 
after  take  a'  constant  in  [5b]  and  [10b].  Combining  [5b]  and 
[lOb[  leads  to  the  result 

1+Ai 

[22] 


wi  -  ^tWt  or  wi/wi0  *  -r  ) 


and,  substituting  [22]  and  [5b]  into  [lib]  leads  to  the  dif¬ 
ferential  relation. 


^dRl  r(l'Al  *  1 


»  R 
dz 


dR 

dz 


[231 


which  has  the  solution, 

Ai 


*  - —  +  K*R  ♦  const.  (Axj^  0) 
Ax 


[2A] 
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and  a'Kiz  =  in  R  +  KaR  +  const.  (A  =  0) 

or ,  (z-z  ) 

a'Ki  — — p-  =  in  R  +  Ka  (ft  -  l)  [24a] 

o 

Substituting  dR/dz  derived  from  [23]  into  [ 5b]  yields  a  relation 
between  ^  and  R, 


ft 

jAi-l) 

R  +  Ka 


[25] 


and,  finally,  it  cay  be  shown  that. 


W  "  _2Ai  Jl+Ai)  (26] 

o  R  +  KaR 

The  type  of  motions  which  ensue  from  this  theory  in  the 
case  of  strong  circulation  are  seen  to  depend  very  much  on  the 
value  of  Ax,  the  constant  appearing  in  the  relatic  for  circula¬ 
tion  change,  and  which  depends  in  part  on  the  way  in  which 
vorticity  is  Ingested  into  the  vortex  pair.  In  fact,  the  asymp¬ 
totic  behaviour  of  the  vortex  pair  changes  radically  as  Ax  varies 
around  the  value  unity.  This  is  demonstrated  in  the  table  below. 
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Asymptotic  Behaviour  (ip  »  l) 


z 

0 

z 

Ax  >  1 

Ai 

R 

Aio'Ki 

-Ax 

Ki  R 

-2Ax 

Kx  R 

-  t* 

Ai  =  1 

n 

Ki/l+Ka 

-2 

Kx  R 

.t* 

Ai  <  1 

R  (Ka/tt'Ki) 

Kx/Ka 

_-(l+Ax) 
Kx/Ka ’R 

t 

Ai  =  0 

R  (Ka/a'Ki) 

Kx/K2 

-1 

Kx/Ka  -R 

.t* 

The  values  of  Ax  >  0,  the  velocity  ratio  is  seen  to  decline, 
so  that  the  strong  circulation  assumption  must  eventually  be¬ 
come  invalid.  The  case  Ai  =  1  yields  results  qualitatively 
similar  to  the  weak  circulation  case.  In  the  case  where  Ax  <  1, 
however,  the  velocity  ratio  increases  the  asymptotic  value  shown 
and,  most  interesting,  the  added  momentum  coefficient  (Kx-Ka^) 
vanishes  asymptotically,  so  that  the  motion  becomes  determined  by 
volume,  vorticity,  and  energy  balances  alone.  Finally,  in  the 
case  where  kx  «  1  (effective  annihilation  of  ingested  vorticity) 
then  asymptotically  the  motion  becomes  determined  by  volume  and 
energy  balances  alone,  yielding  z  ~  t2' 
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III.  COMPARISON  WITH  EXPERIMENT  (HOMOGENEOUS  MEDIA) 

In  Figures  5  and  6  are  shown  data  from  actual  experiments 
on  two-dimensional  vortex  pair  motions  in  homogeneous  fluids. 

Most  of  the  data  shown  were  obtained  in  experiments  carried  out 
in  our  own  laboratory.  Suffice  it  here  to  show  a  schematic  of 
the  facility  which  was  used.  Fig.  7,  and  to  show  Table  1,  in 
which  the  properties  and  characteristics  of  the  experimental 
vortex  pairs  are  listed. 

Most  significant,  we  found  in  our  experiments  and  from  the 

data  of  Richards  (1965)  that  the  measured  variation  of  vertical 

velocity  and  pair  radius  (two-dimensions)  conformed  more  closely 

-1  i 

to  the  law  W  ~  R  or  z  ~  t  than  to  the  law  derived  in  the  past 

by  others  and  which  is  based  on  complete  similarity  and  momentum 

2  \ 

conservation]  i.e.  W  «•  R  or  z  ^  t  .  A  test  of  the  simple  con¬ 
servation  of  momentum,  W  ~  R~2,  using  a  typical  trajectory  is 
illustrated  In  Fig.  5  and,  similarly,  in  Fig.  6  it  Is  shown  that 
the  trajectories,  so  far  as  they  have  been  observed  experimentally, 
conform  more  closely  to  the  asymptotic  law  derived  earlier  for 
the  case  of  strong  circulation,  utilizing  a  small  value  of  Ax 
(0  <  Ax  <  2). 

In  the  case  of  strong  circulation,  the  radius  grows  in  a 
linear  fashion  asymptotically,  but  the  theory  predicts  that 
during  the  initial  phases  of  the  motion  the  quantity  dR/dz  is 
less  than  its  asymptotic  value.  A  similar  behaviour  was  observed 
In  our  experiments,  see  Fig.  4.  The  matching-up  of  these  ob¬ 
served  trajectories  with  the  theory  offers  an  opportunity  to  de¬ 
termine  some  of  the  constants  of  the  theory.  For  this  purpose 
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we  assume  to  begin  with  that  Ai  =  0,  since  the  comparison  be¬ 
tween  observed  end  theoretical  trajectories  suggests  a  small 
value.  In  this  case, 

z-z 

(a'Ki)  =  In  R  +  Ka  [R-l]  [24a] 

o 

The  trajectory  according  to  Eq.  [24a]  with  Ks  =  1/2  and  a'Ki  =  .38 
is  shown  in  Fig.  4.  A  fair  fit  with  the  experiments  has  been 
achieved,  and  noticeably  better  than  is  possible  with  any  linear 
trajectory. 

The  strong  circulation  theory  thus  explains  the  two  impor¬ 
tant  features  of  vortex  pair  behaviour  which  cannot  be  explained 
by  the  usual  theory  of  complete  similarity.  These  features  are: 

(i)  t’-’e  tendency  for  forced  vortex  trajectories  (homogeneous 

1  A 

flow)  to  more  nearly  follow  the  law  z  ~  t?  rather  than  z  ~  ts, 
and  (ii)  the  tendency  for  the  entrainment  coefficient  (dR/dz)  to 
grow  during  the  initial  phase  of  the  motion.  These  results  sug¬ 
gest  that  in  vortex  motions  in  homogeneous  flows  the  internal 
velocity  scale  grows  steadily  relative  to  the  translational 
(external)  velocity,  the  ratio  approaching  a  value  considerably 
larger  than  unity,  while  at  the  same  time,  the  virtual  momentum 
coefficient  associated  with  the  vortex  motion  approaches  the 
value  zero.  The  data  also  suggest  that  vorticity  ingested  from 
around  one  half  of  the  vortex  pair  is  almost  annihilated  through 
mixing  with  vorticity  ingested  from  the  opposite  side. 
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IV.  SIMPLIFIED  THEORY  {VORTEX  PAIR  MOTION  IN  STRATIFIED  MEDIA) 


Convected  masses  In  nature  are  often  rising  or  falling  in  a 
medium  of  varying  density,  as  in  the  case  of  a  chimney  plume  pro¬ 
jected  upwards  into  a  stable  atmosphere.  The  latter  may  be  char¬ 
acterized  by  a  characteristic  time  (the  Vaissala  period),  1/ 

1  dpe 

- - tt  and  p  is  the  potential  density  of  the  atmo- 

»  a 


where  a  = 


sphere.  The  same  definition  can  be  used  to  characterize  any 
density  stratified  media. 


The  motion  of  the  convecting  mass  may  also  be  characterized 
at  any  instant  by  the  time,  R/W.  It  is  almost  apparent  that  when 
the  latter  time  is  long  in  comparison  to  the  Vaissala  period  that 
the  effect  of  stratification  will  dominate,  and  conversely.  That 
is. 


Effect  of 

Stratification 

Vanishes 


decreasing 

<3 - 


increasing 

- o 


Stratification 

Dominates 


Quite  clearly,  too,  as  the  motion  proceeds  in  time,  the  ratio  R/W 
increases  continuously,  so  that  stratification  must  eventually 
dominate.  When  this  happens,  the  vertical  motion  of  vortex-pairs 
may  become  oscillatory,  and  is  accompanied  by  the  collapse  and 
horizontal  spreading  of  the  convected  mass,  as  illustrated  in  Fig¬ 
ure  8.  This  behavior  is,  of  course,  not  consistent  with  similarity 
either  complete  or  of  the  kind  assumed  in  the  preceding  section. 


It  is  sometimes  desirable  .  j  be  able  to  estimate  the  tra¬ 
jectory  of  a  vortex  pair  while  it  is  rising  in  a  stratified  media 
and  particularly  to  predict  the  maximum  height  of  rise  and  the  time 
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required  to  reach  the  maximum.  For  this  purpose,  we  adopt  here 
a  simplified  theory  based  essentially  on  the  assumption  of  strong 
circulation  and  annihilation  of  ingested  vorticity.  In  fact,  the 
particular  assumptions  adopted  would  apply  if  the  velocity  ratio, 
had  already  closely  approached  its  limiting  value.  These  as¬ 
sumptions  are:  (i)  the  motion  is  determined  by  conservation  of 
volume,  mass,  and  energy  (neglecting  vorticity  and  momentum); 

(ii)  complete  similarity  (dR/dz  =  0,  a  constant).  For. further 
justification  of  these  assumptions  we  shall  depend  finally  upon 
a  comparison  between  theoretical  predictions  and  the  results  of 
systematic  experiments. 

The  energy  balance  is  expressed  as  follows. 


bE  _  j*. 

dt  8t 


|  (u2+v2+w2 )  +  (p-pe)  gC  d£dqd£  = 


See  Figure  14  for  nomenclature. 


(Rate  of  Dis¬ 
sipation  of 
Energy) 

[27] 


For  a  self-similar,  self-convecting  flow,  the  dissipation 
of  kinetic  energy  per  unit  volume  which  occurs  due  to  the  action 
of  turbulent  shear  stresses  must  for  dimensional  reasons  be  of 
the  form. 

Dissipation  pW3  r 1 

Unit  Volume  *  R 


As  a  result,  the  energy  balance.  Equation  [27],  for  a  self- 
convecting  mass  in  a  homogeneous  medium  of  the  same  density  takes 
the  form. 


K  d(W2R2+J)  „  W3  „2+J 

2  dt  =  ’  °D  T  R 


0 


[29] 
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where 

J  =  0  (two-dimensions) 
j  =  1  (axisymmetric ) 

and  K  is  a  constant  of  virtual  energy,  defined  by  the  identity. 


ug+  v2+  W2 
2 


d^dr)dC  =  |  W2R2+J 


[30] 


and  where  Cp  is  a  dissipation  coefficient. 

Making  use  of  [29],  together  with  the  relationship  R  =  £z, 
it  may  be  shown  that  the  height  of  rise  follows  the  law. 


z  (2+J/2  +  V0K)  ^  t 


[31] 


in  a  medium  of  uniform  density  with  Ap  =  0.  The  dissipation 
coefficient  in  nature,  D  =  may  be  determined  by  a  com¬ 

parison  between  theoretical  trajectories  such  as  given  by  [31] 
and  observations  of  vortex  pair  rise  in  homogeneous  media.  As 
shown  in  Fig.  6,  such  a  comparison  leads  to  the  conclusion  that 
D  is  quite  small  (D  <  0.2). 

The  trajectory  given  by  [31]  may  be  compared  to  the  law 
which  would  apply  if  momentum  were  conserved, 

z3+J  «  t  [32] 

which  coincides  with  [31]  only  if  D  ■  1.  The  variance  of  ob¬ 
served  trajectories  from  the  momentum  law  [32]  is  clearly  seen 
in  Figs.  5  and  6. 
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Volume  conservation  in  a  self-similar  flow  leads  to  a 
linear  relation  between  the  nominal  radius  of  the  mass  and  the 
height  of  rise  from  the  virtual  origin  z  =  0: 

R  =  jSz  [33] 

Conservation  of  mass  takes  the  form 


_d_ 

dt 


4 

3 


ttR 


2+j 


Pi 


2ttR1+J  peWiS 


[3k] 


where  pg,  is  the  density  of  the  surrounding  fluid  at  any  given 
height  z  and  p  is  the  average  density  within  the  rising  mass, 
defined  by 


]  d 

-  ^R2+J  [ p± (z )  -  pe(z)]  =  /  (p-pe)  d^(dT!)JdC  ,  [35] 


where  the  integration  is  taken  over  the  entire  volume  of  the 
rising  mass. 

The  formulation  of  conservation  of  energy  is  based  upon 
[27]  and  [28] 


d_ 

dt 


K 

2  PiW 


*r2+j 


+  k(pi-  pe)gzR2+J 


wf_ 

D^i  R 


R2^ 


[36] 


where  W  and  R  are  the  observed  and  measured  gross  properties  of 
the  rising  mass  while  K  and  k  are  the  coefficients  of  the  virtual 
kinetic  and  potential  energies,  respectively,  defined  in  two  di¬ 
mensions,  e.g.,  by 
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KPi  -  =  /  §  (ua+  w2)  d£d£  [3Tl 

and 

k(p1-  pg)  zR2  =  /  (p-pe)(z+0  d£d£  [38] 

In  most  practical  Instances  where  one  is  dealing  with  a 
mass  of  fluid  convected  through  a  homogeneous  or  stratified 
medium  such  as  the  ocean  or  the  atmosphere*  the  difference  be¬ 
tween  the  densities  of  the  convected  and  surrounding  masses  is 
very  small;  that  is  Ap/pg  «  1*  being  usually  of  the  order  of 
10"3,  and  therefore  p^/pe  can  be  taken  as  1.  This  assumption, 
frequently  referred  to  as  the  Boussinesq  approximation,  see 
Phillips  (1966),  will  be  used  throughout  the  analysis  presented 
herein. 

Using  the  Boussinesq  approximation  and  the  identity  dz  =  wdt, 
the  three  conservation  statements,  Equations  [33l*  [34]  and  [36], 
may  be  reduced  to  the  following  form  in  the  case  of  a  planar 
motion: 

R  =  ]3z  [39] 


+  az 


z  ss  0 


[41] 


1 
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where  a  =  -  - - —  and  Ap  =  (p  -  pg). 

Finally,  explicit  general  solutions  of  [40]  and  [4l]  may 
be  found.  They  are: 


Afi.  = 

|A£.|  .  !f° 

^e 

L  '0  3  J 

az 

--2  ,  _ o  - 

z  +  -y  z 


w 

w 


1  - 


+  B 


(1  +  2D)  1  +  D/2  1  +  2D 


-  -n 
z 


r  (a  +  b)  i  -_i  r _ b__J 

[(1  +  2D)  J  2  '  [(1  +  D/JJ 


where 


A  = 


2k 

K 


G  and  G  s 


2k 

3K 


Bs  ”  S  and  S  = 


isi 

W  4  f>e 


z  2ag 
o 


W 


D  = 

^  ■“  Atr 


D 

0K 


and 


n  s  2 ( 1+D ) 

z  s  z/z 


[42] 


[431 


[44] 

[45] 

[46] 

[47] 

[48] 
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The  parameter  G  is  a  measure  of  the  initial  buoyancy  of  the 

converted  mass  relative  to  its  initial  momentum.  G  is  taken  as 

positive  when  a  net  buoyancy  force  is  acting  on  the  mass,  i.e. 

Ap  <  0.  S  is  a  measure  of  the  added  buoyancy  which  would  result 

from  moving  the  convected  mass  a  vertical  distance  zq  through  a 

stratified  medium.  Since  Vli  is  the  frequency  with  which  a 

finite  volume  of  fluid  of  given  density  would  oscillate  in  a 

stratified  medium,  often  referred  to  as  the  Vaisala  frequency, 

the  parameter  S  can  also  be  considered  as  the  square  of  the  ratio 

of  the  characteristic  time  of  the  convected  motion,  z  /W  ,  and 

o  o 

the  reciprocal  of  the  Vaisala  frequency. 


The  maximum  height  of  the  rising  mass,  reached  at  the  point 
where  W  =  0,  is  according  to  [43],  given  by  the  solution  of  the 
following: 


2+n 


r  b  i 

Z  \ 

max  1 

r 

A  +  B  1  1 

z  ) 
max  j 

[(l  +' D/2)J 

z°  ) 

i 

1  +  2D  ]  | 

0  J 

+  B 


1  +  D/2  1  +  2D 


A 

(1  +  2D) 


[49) 


It  is  of  Interest  to  consider  certain  special  cases: 

1.  A  mass  rising  in  a  homogeneous  medium  with  the 
same  density  as  itself;  i.e.,  A  -  0  and  B  -  0.  Then, 


z_ 

z 


-n/2 


(90] 


or, 


l+n/2 


z 

z 


1  + 


1  +  11 
2 


W  t 
o 


(51) 
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This  result  suggests  how  to  estimate  the  dimensionless  quantity 
n  (or  D)  through  the  analysis  of  the  trajectories  of  rising  masses 
in  this  special  case. 

2.  A  mass  with  initial  density  difference  rising  in 
a  homogeneous  medium;  i.e.,  B  =  0.  Then, 


If  A  >  0,  then  no  maximum  height  is  reached,  but  if  A  <  0, 

1 

(1+2D) 

1  +  2D  +  |a| 

- M -  [53] 

1*1  J 

The  predicted  rise  of  the  mass  as  a  function  of  time  and  the 
maximum  rise  of  the  mass  for  a  range  of  values  of  A(<0)  and  D,  as 
obtained  from  Equations  [52]  and  [53],  are  presented  in  Figures 
9  and  10.  These  figures  demonstrate  clearly  the  effect  of  the 
(negative)  initial  buoyancy  and  energy  dissipation  parameters  on 
the  time  history  of  an  impulsively  started  rising  mass  moving 
through  a  uniform  surrounding  fluid  of  smaller  density. 

2a.  The  same  case  as  above  but  for  W  -  0  and  A  >  0. 

o 

First  of  all,  [43]  may  be  rewritten: 


: 
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or  in  this  case 


W2.  2  =  _ i2kg _  /Afi. 

02 (1  +  2D)K_  '  pe 


3.  A  mass  with  no  initial  buoyancy  rising  in  a  strati¬ 
fied  medium,  i.e.,  A  =  0.  Then, 


JL  n  .  B  1  1  -  -n  B  -_i 

Wq  |  1  +  D/2  "  1  +  2D  |  Z  (i  +  D/2)  ‘  Z 


(1  +  D/2. 


i2  [551 


and  the  maximum  rise  of  the  mass,  as  a  function  of  B,  is  obtained 
from  Equation  [55]  by  setting  W  =  0. 

Approximate  integrated  solution  for  the  height  of  the  con- 
vected  mass,  z,  as  a  function  of  time  can  be  readily  obtained 
from  Equation  [55]  whenever  D  <<  1  and  B  «  ’ .  When  these  two 
requirements  are  satisfied,  Equation  [ 55 i  can  be  rewritten  as 


*  z~ 2  -  Bz  2 


[56] 


and  upon  integration  we  obtain 


i2  2  *  S 


in  Jsin-1  V"b  +  2  VSj 


where  we  have  normalized  the  time  t  according  to 


(57) 


[58] 
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Equation  [57]  is  particularly  useful  for  the  approximate 
determination  of  the  maximum  rise  of  a  mass  convected  in  a  strati¬ 
fied  medium  and  the  time  at  which  this  maximum  height  is  reached, 

t  .  For  the  maximum  height  we  find 
max 


z 

max 


[59] 


and  the  time  required  to  reach  this  height  is  given  by 


t 

max 


7T  1 

4  Vi'  ? 


[60] 


This  last  result  is  especially  interesting.  In  experiments 

on  convected  masses  of  fluid  moving  through  a  density  stratified 

medium  it  was  frequently  observed  that  the  time  it  takes  the 

mass  to  reach  its  maximum  height  is  inversely  proportional  to 

the  Vaisala  frequency  of  the  stratified  medium,  i.e., 

t  Vag  =  const.  Equation  [60]  is  Just  a  statement  of  this  same 
max  v 

fact  for  small  values  of  B  (as  usually  exist  in  nature),  since 

tmaxVi  =  Vii/ii  (tmaxV^i). 

V.  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL 
RESULTS  (STRATIFIED  MEDIA) 

In  our  experimental  investigation  we  have  studied  the  motion 
of  impulsively  started  rising  masses  (or  vortex-pairs)  in  both 
uniform  and  density-stratified  surroundings.  According  to  the 
theoretical  considerations  presented  in  the  previous  section, 
the  time-histories  of  these  motions,  expressed  in  appropriate 
non-dimensional  terms,  are  determined  by  three  parameters,  A,  B 
and  D ;  i.e.,  for  given  values  of  these  parameters  the  rise  and 
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growth  of  the  convected  mass  as  a  function  of  time  can  be  pre¬ 
dicted.  A  and  B  are  determined  by  G  and  S  and  by  a  third  param¬ 
eter  which  is  the  ratio  of  the  virtual  kinetic  and  potential 
energy  coefficients,  K/k,  according  to  Equations  [44]  and  [45]. 

While  the  parameters  G  and  S  are  determined  in  each  case  by 
the  initial  conditions  of  the  rising  vortex-pair,  there  is  no 
practical  way  for  determining  a  priori  the  values  of  D  and  K/k. 
These  latter  parameters  can  be  determined  only  by  comparing  cer¬ 
tain  sets  of  experimental  results  with  corresponding  theory. 

A  series  of  experiments  on  the  motion  of  vortex-pairs  in 
a  homogeneous  medium  of  the  same  density.  Series  III  (see  Table  l), 
where  the  parameters  G  and  3  (and  therefore  also  A  and  B)  are 
identically  equal  zero,  may  be  used  for  determining  the  dissipa¬ 
tion  parameter  D.  The  rise  of  the  vortex-pairs  in  this  case  is 
predicted  by  Equation  [51]  and  is  graphically  depicted  in  Fig¬ 
ure  9  (with  A  =  0).  The  actual  predicted  rise  of  the  convected 
mass  depends  on  the  numerical  value  of  the  parameter  D  (or  n). 

In  Figure  6  are  shown  a  comparison  between  experimental  and 
theoretical  results  on  the  rise  of  impulsively  started  masses  in 
a  uniform  medium  of  the  same  density.  In  a  log-log  plot,  the 
slope  of  the  trajectory  for  large  values  of  (WQt/zo)  should  be 
equal,  according  to  our  analysis,  to  1/2  +  D,  and  it  can  be  used 
therefore  for  determining  the  value  of  the  dissipation  parameter 
D  associated  with  the  motion  of  the  rising  mass.  Included  in 
Figure  6  are  the  experimental  results  of  Richards  ( 1 965 ) j  on  the 
rise  of  two-dimensional  puffs  in  homogeneous  surroundings.  The 
best  agreement  with  all  experimental  results  is  obtained  when 
we  choose  D  =  0.2. 
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The  numerical  value  of  K/k  enters  into  the  analysis  only 
when  there  is  an  initial  difference  between  the  rising  and  sur¬ 
rounding  fluid  densities  or  when  the  surrounding  fluid  is  strati¬ 
fied.  This  value  will  be  also  determined  from  a  comparison  of 
some  experimental  and  theoretical  results.  For  a  vortex-pair 
convected  in  a  density-stratified  medium  we  found  earlier  that, 
for  sufficiently  small  values  of  the  parameters  A,  B  and  D,  the 
time  it  takes  the  mass  to  reach  its  maximum  height  is  inversely 
proportional  to  the  Vaisala  frequency  and  is  given  by 

“Vi?  [6l! 

This  value  decreases  only  very  gradually  as  the  value  of  B  (or 
S)  sc  chat  Eq.  [6l]  is  very  useful  for  the  experimental  deter¬ 
mination  of  K/k. 

In  Figure  11  the  value  of  the  product  (t  o  V&g)*  as  mea- 

max 

sured  in  the  experiments  of  Series  I  and  II,  is  presented  as  a 
function  of  the  stratification  parameters  S;  the  initial  condi¬ 
tions  for  each  experiment  presented  in  the  Figure  are  included 
in  Table  1.  There  are  certain  inherent  inaccuracies  in  the  ex¬ 
perimental  determination  of  tmgx  which  explain  the  scatter.  Also 
shown  in  Figure  11  are  the  asymptotic  solution  for  the  maximum 
rise  time,  Equation  [60],  and  the  exact  solution,  according  to 
Equation  [43],  with  D  =  0.2  and  G/S  =  -0.715-  Close  agreement 
between  theoretical  and  experimental  results  was  obtained  when 
we  used  3K/2K  =  6  or  K/k  =4.  We  have  included  in  the  same  fig¬ 
ure  the  experimental  results  from  Series  III  which  had  a  markedly 
different  value  of  j3;  these  too  were  found  to  agree  very  closely 
with  the  theoretical  prediction  based  on  D  =  0.2  and  K/k  =  4, 
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indicating  that  the  dependence  of  these  latter  two  parameters  on 
£  is  probably  very  weak.  We  have  used  these  values  of  D  and 
K/k  for  all  subsequent  comparisons  of  experimental  and  theoreti¬ 
cal  results.  The  coefficient  of  virtual  potential  energy  cannot 
be  much  different  from  unity,  according  to  its  definition  in 
Equation  [38].  The  total  kinetic  energy  of  a  rising  vortex-pair 
was  thus  found  to  be  about  four  times  larger  than  the  kinetic 
energy  associated  with  its  linear  convection  alone,  an  indication 
of  the  intensity  of  motion  inside  the  vortex-pair,  which  con¬ 
tributes  to  its  total  kinetic  energy.  This  finding  lends  im¬ 
portant  support  to  the  strong  circulation  assumption. 

A  comparison  between  the  predicted  and  actual  maximum  heights 
of  rise  of  a  vortex-pair  convected  in  a  linearly  density-strati¬ 
fied  medium  is  shown  as  Figure  12.  The  figure  includes  a  pre¬ 
diction  based  on  the  simplified  asymptotic  solution  for  small  B 
(or  S),  Equation  [59],  and  a  prediction  obtained  from  the  exact 
solution  of  Equation  [49]  for  zm  .  Generally  there  is  good 
agreement  between  the  experimentally  measured  maximum  height  of 
vortex-pairs  and  the  exact  theoretical  solution. 

Finally,  in  Figure  13,  we  compare  the  measured  trajectories 
(height  versus  time)  of  vortex-pairs  with  their  theoretically 
predicted  trajectories.  The  vortex-pairs  included  in  the  Fig¬ 
ure  all  had  different  starting  conditions  and  they  were  moving 
through  media  with  different  density-stratifications.  However, 
their  trajectories,  as  depicted  in  the  figure  are  shown  to  de¬ 
pend  only  on  the  values  of  the  two  lumped  parameters  G  and  S 
which  combine  their  starting  conditions  with  the  properties  of 
the  surrounding  medium. 


HYDRONAUTICS,  Incorporated 

-30- 

The  fact  that  trajectories  of  different  vortex-pairs  are 
grouped  according  to  their  G  and  S  values  confirms  the  validity 
of  the  scaling  laws  and  scaling  parameters  used  herein,  while 
the  agreement  obtained  between  the  experimental  and  theoretical 
trajectories  lends  further  support  to  the  validity  of  the  simpli¬ 
fied  theory  presented  here  for  the  motion  of  vortex-pairs  in 
stratified  media. 


VI.  SUMMARY  AND  CONCLUSIONS 

A  theory  for  the  motion  of  two-dimensional  turbulent  vortex- 
pairs  in  homogeneous  media  has  been  developed  based  on  separate 
velocity  scaling  of  the  internal  and  external  flow  fields  in¬ 
volved  in  the  motion.  These  two  flow  fields  are  depicted  to  be 
separated  by  a  thin  region  of  high  shear,  which  also  forms  the 
boundary  of  the  captured  mass.  The  theory  takes  into  considera¬ 
tion  variations  in  volume,  circulation,  momentum,  and  energy  in 
the  flow  field.  The  ratio  of  internal  to  external  velocity 
scales,  ip,  is  introduced  as  an  important  variable.  The  virtual 
momentum  coefficient  is  shown  to  be  linear  in  ip,  of  the  form 
Ki-K  zip. 

The  theory  is  specifically  derived  for  the  two  limiting 
cases  of  weak  and  strong  circulation.  In  the  former  case, 
ip  -*•  1  and  the  entrainment  is  weak;  the  asymptotic  behavior  of 
the  trajectory  is  z  ~  t5  Just  as  predicted  by  the  usual  theory 
based  on  complete  similarity  and  momentum  conservation. 

In  the  case  of  strong  circulation,  ip  »  1,  the  asymptotic 
behavior  of  the  trajectory  depends  very  much  on  the  way  in  which 
vorticity  from  the  shear  layer  is  ingested  into  the  vortex  pair. 
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/ 

/ 

In  the  case  where  the  shear  ia'yer  from  opposite  sides  is  in¬ 
gested  in  such  a  way  as  to  cause  annihilation  of  the  ingested 

l. 

vorticity,  then  the  asymptotic  trajectory  is  z  ~  t2.  Under  the 
same  conditions,  the  velocity  ratio,  f,  increases  toward  the 
asymptotic  value  K1/K2  so  that  the  virtual  momentum  coefficient 
tends  to  zero.  As  a  result,  the  asymptotic  motion  assuming 
vorticity  annihilation  corresponds  to  a  motion  with  complete 
similarity  and  with  energy  conservation.  The  ratio  of  growth 
of  the  pair  radius  with  height  is  shown  to  increase,  approaching 
a  linear  relation  asymptotically. 

Systematic  experiments  have  been  carried  out,  and  the  re¬ 
sults  for  rise  versus  time  and  radius  versus  height  are  compared 
with  the  theory.  They  lend  strong  support  to  the  strong  circu¬ 
lation  theory  and  further  suggest  that  ingested  vorticity  is  to 
a  large  degree  annihilated. 

Based  on  these  findings  for  the  case  of  homogeneous  flows, 
a  simplified  theory  is  derived  for  the  rising  motion  of  vortex 
pairs  in  stratified  media.  The  assumptions  of  the  theory  are: 

(i)  the  motion  is  determined  by  conservation  of  volume,  mass, 
and  energy  (neglecting  vorticity  ar.d  momentur.);  (ii)  complete 
similarity  (dR/dz  =  0,  a  constant).  General  laws  of  motion  in 
stratified  media  have  been  derived  and  solutions  given;  particu¬ 
larly  interesting  cases  are  discussed  ir.  detail. 

Motions  in  stratified  media  were  shown  to  depend  on  four 
non-dimensional  parameters.  Two  of  these  depend  upon  the  initial 
conditions  of  the  motion  and  the  stratification  of  the  media. 

The  other  two  are  Inherent  in  the  details  of  the  motion  and  had 
to  be  determined  from  experiments;  one  of  tnese,  the  dissipation 
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parameter  D  =  Cjy//SK  was  found  to  be  0.2  while  the  other,  the 
ratio  of  virtual  kinetic  and  potential  energy  coefficients  K/k 
was  found  to  be  4.  On  the  basis  of  these  numbers  it  may  be 
concluded  that  the  dissipation  rate  is  small  and  that  the  con- 
tribu" ion  of  internal  motions  to  the  overall  kinetic  energy  is 
large. 

The  experiments  confirmed  the  environmental  scaling  param¬ 
eters,  which  were  used  to  collapse  data  taken  under  differing 
conditions.  Good  agreement  was  found  between  predicted  and  ob¬ 
served  trajectories.  Particularly  good  agreement  was  found  for 
the  maximum  height  of  rise.  The  time  required  to  reach  maximum 
height  was  found  to  be  inversely  proportional  to  the  Vaisala 
frequency,  'J'ag,  and  was  given  approximately  by  t  Vag  =  1.8, 
in  good  agreement  with  the  theory.  In  general  the  experiments 
confirmed  the  utility  of  the  simplified  theory  for  predictions 
of  the  motion  of  vortex  pairs  in  stratified  media.  This  theory 
has  been  utilized  elsewhere  for  the  prediction  of  the  behavior 
of  chimney  plumes  rising  into  a  stable  atmosphere,  with  very 
good  agreement  between  the  theory  and  full  scale  observations, 
Tulin  and  Shwartz  (1970),  and  also  with  excellent  correlation 
with  experiments  to  the  penetration  of  a  density  discontinuity 
by  a  turbulent  vortex-pair,  Birkhead,  Shwartz,  and  Tulin  (1969). 
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FIGURE  I  -  EXAMPLES  OF  VORTEX- PAIRS 
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-  THE  VERTICAL  VELOCITY  VS.  RADIUS  OF  A  VORTEX -PAIR  MOVING 
IN  A  HOMOGENEOUS  FLUID 

^7 


HYDRONAUTICS,  INCORPORATED 


TIME,  Wot/z0 


FIGURE  6  -  THE  RISE  OF  VORTEX -PAIRS  IN  A  HOMOGENEOUS  MEDIUM;  COMPARISON 
OF  EXPERIMENT  AND  THEORY 
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FIGURE  7  -  EXPERIMENTAL  FACILITY  FOR  STUDYING  VORTEX  PAIR  MC 
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FIGURE  8  -  TYPICAL  TRAJECTORY  OF  A  VORTEX-PAIR  IN  A  LINEARLY 
DENSITY  STRATIFIED  MEDIUM 


FIGURE  9  -  THE  RISE  OF  AN  IMPULSIVELY  STARTED  HEAVY  MASS  OF  FLUID 
IN  UNIFORM  SURROUNDINGS  (B  =  0),  THEORY 
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FIGURE  11  -  TIME  OF  MAXIMUM  RISE  FOR  VORTEX-PAIRS  CONVECTED  IN  A  DENSITY-STRATIFIED 
MEDIUM;  COMPARISON  OF  EXPERIMENT  AND  THEORY 
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(1)  Eq.  21,  p.  12  should  be  renamed  Eq.  21a. 

In  this  equation,  in  the  left-hand  side,  a  should 
be  replaced  by  a',  and  a  (’Jr)  in  the  preceeding 
sentence  should  be  replaced  by  a'. 

(2)  Eq.  23,  p.  12,  Right-hand  side  should  be  multiplied 

by  u,  where  q  =  (M  / p )/ (W  R*),  i.e.,  Eq.  23  should  now 

Z  X  o 

read 


K, 


Kj.  dR  )  _  U  dR 


ds 


a 1 

dz 


C23] 


(3)  Eq.  24,  p„  12,  should  be  replaced  by 


a  1 KT  z 
U 


+  const. 


(V  °) 


[24] 


And  the  following  equation,  on  p.  13,  by 


a'K,  z  —  .  KaR  . 

- 1 —  =  in  R  +  — 8 —  +  const. 

u  u 


(A  =  0) 
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(8)  Eq.  25a*  p.  16,  should  be  replaced  by 


.  T r  Z-Z  IT  — 

SL£l  =  _£  =  m  R  +  £*-  [R  -  1] 

u  R0  u 


[25a] 


